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Tuning of the Gap in a Laughlin-Bychkov-Rashba Incompressible Liquid 
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We report on our investigation of the influence of Bychkov-Rashba spin-orbit interaction (SOI) on 
the incompressible Laughlin state. We find that experimentally obtainable values of the spin-orbit 
coupling strength can induce as much as a 25% increase in the quasiparticle-quasihole gap E g at 
low magnetic fields in InAs, thereby increasing the stability of the liquid state. The SOI-modulated 
enhancement of E g is also significant for v — 1/5 and 1/7, where the FQH state is usually weak. 
This raises the intriguing possibility of tuning, via the SO coupling strength, the liquid to solid 
transition to much lower densities. 

PACS numbers: 73.43.-f,73.21.-b 



The ground state of a two-dimensional electron gas 
(2DEG) in the presence of a high perpendicular mag- 
netic field is well known to appear as a multitude of 
highly correlated incompressible fractional quantum Hall 
effect (FQHE) states at a few special values of the Landau 
level filling factors Q, • It is also well established that 
the strongest effect appears at the lowest Landau level 



filling factor 



i) that has the largest quasiparticle- 



quasihole gap 0,12- The gaps are much smaller for v < 
and as a consequence, nature of the electron states in 
that regime has remained a challenge until now. This is 
due to the fact that in the low density regime, the in- 
compressible liquid is expected to undergo a phase tran- 
sition to a crystalline state Q- However, a definitive 
conclusion about the onset of this quantum phase tran- 
sition has remained elusive because experimentally one 
observes weak effects for filling factors \ etc., and the- 
oretically, one compares two very small energies (ground 
state) in order to determine which phase is energetically 
favored. For more than two decades, investigations of the 
FQHE have focused largely on 2DEGs that are embedded 
in GaAs heterostructures, where spin-related effects are 
small (though important 4j) compared to other effects, 
because of the small value of the ^-factor of electrons in 
GaAs |5(. However, studies of the spin-orbit (SO) cou- 
pling in a 2DEG within an InAs (or InSb) quantum well 
with very large g values, are at the cusp of a rapid ad- 
vance, due largely to their relevance to spin transport in 
low-dimensional electron channels ;6:j . In order to investi- 
gate the influence of SO interaction on the incompressible 
Laughlin states, we have carried out the well established 
finite-size studies in a periodic rectangular geometry, but 
for the first time with the SO coupling included in the 
Hamiltonian. We find that as the SO coupling strength 
is increased, there is an increase in the quasiparticle- 
quasihole gap, indicating that the incompressible state 
can be rendered more stable by appropriate tuning of 
the SO coupling. This might prove to be particularly 
useful in the low-electron density regime where, as de- 
scribed above, the FQHE state is usually very weak for 
conventional systems. 



For a 2DEG in the xy plane with a magnetic field B 
along the z direction, in the Landau gauge [with the 
choice of vector potential A = (0, Sac, 0)], the single- 
particle states and the corresponding energies are ob- 
tained by solving the one-electron Hamiltonian 



H 



(p + eA) 2 
2m* 



'-[a x (p + eA)] z +gn B Ba z (1) 



that includes the Bychkov-Rashba term and the Zee- 
man term. Here p is the momentum operator, a is the 
SO coupling strength and cr = (a x , o v ,a z ) are the Pauli 
spin matrices. Experimental values of the SO coupling 
constant lie in the range of 5-45 meV-nm [8|. The high 
values of a are deduced from magnetotransport experi- 
ments where the SO interaction is tuned for a fixed carrier 
density in a 2DEG by using gate electrodes in a square 
asymmetric InAs quantum well Q. We will therefore fo- 
cus our investigation on InAs, as it represents the most 
promising material, as yet, for achieving large SO cou- 
pling. 

We solve the Schrodinger equation 



Hip = Eip 



(2) 



in a rectangular geometry with supercell sides L x and 
L y (i.e., aspect ratio A = L x /L y ) and expand the single- 
particle wavefunctions ipk y (r) as a superposition of solu- 
tions of the Hamiltonian in the absence of SO interaction 
i 



(p n (x - X )C° |(T>/a 



(3) 



where 



4> n (x - ar ) = n e- ix - Xo)2/2l °H n [(x - x )/l ]/^M> 

is the usual solution to the harmonic oscillator problem, 
with H n (x) the Hermite polynomial, Iq = {h/m*uj c ) 1 / 2 
the radius of the cyclotron orbit with frequency lu c = 
eB/m* and center xq — k y lg, n is the Landau level index, 



2 



A, = l/y/2 n n\, and 



a = up, dn 



' 1 



is the electron spinor. 

Substituting Eq. into the Schrodinger equation > 
multiplying both sides by <pi(x — Xq) and integrating over 
x, we obtain a system of equations 



%{a/l Q )y/%C^ l + [{I + l/2)fiw c + E d ]Cf n = 

[(Z + l/2)hw c + E U ]C? P - t(a/l )^2(l + l)C^ 1 = 

1 = 0,1,2,... (4) 

whose solution yields @ 

(l/2fej c + £ d )Cf 1 = 0, s = 0, 



(s- l/2)/k<; c + .E„ -i(a/l )V2s 
i(a/l )V2s (s+l/2)twj c + E d 



= 0, 



s = 1,2,3,..., where -E u = .g/is-B — £ and = 
—gfisB — E. Corresponding to s = there is only one 
level, the same as the lowest Landau level without SO 
interaction, with energy 

E = l/2hw c - g\x B B 

and wavefunction 

4>o,k y = e tkyV Mx ~ x o)/y% 



For all other values of s ^ there are two branches of 
levels @ 

, J+ ^1 (-iDs<t>s-i{x-x Q )\ () 




FIG. 1: Ground state energy per particle as a function of fill- 
ing factor v for 4 different values of the SO coupling constant 
a — 0, 10, 20, 40 calculated for 4 electrons in the lowest two 
Landau levels in InAs (B = 10 Lesla). 



and 

^S,ky 

with energies 



(f) s -l(x - Xq) 

-iD s 4> s {x - x ) 



E: 



skv r ±\/E 2 + 2sa 2 /l 2 . 



Here 



2sa/l 



E + ^E 2 + 2sa 2 /l 2 



(6) 



(7) 



(8) 



and A s = 1 + D 2 . From Eq. JSJ and |JB} we see that 
SO interaction couples two Landau levels. While previ- 
ous works were restricted to the study of the single- 
particle states, we use these equations as a starting point 
for our exact many-body treatment. Applying periodic 
boundary conditions, we obtain (k y — xq/Iq) 

x a = Xj = 2i:l 2 j/Ly, L x = 2-Kl\mjLy, 

and consequently, 



1 



yj y/Trl LyA, 
(Xj + nL x — x) 2 



^exp 



i(Xj + nL x )- 2 



r 2 




V 7 y/nloLyA, 
(Xj + nL x — x) 2 



i(Xj + nL x )-^ 



2/ 2 



(Xj+nL x —x) 



Ps-lHs-l ( 

- iDs fi s H s ( v ' ) 



We then build the antisymmetrized products (Slater de- 
terminants) using ip + and ip~ as a complete basis for the 
many-body wavefunction expansion 

* = ^P(i l! i 2 ,...,i n )at i at 3 ...atj0) 

where i k = (sk,jk,Vk), &k = ± and V(h, 12, ■ ■ ■ , in) is 
the antisymmetrization operator. 

The many-body Schrodinger equation was then solved 
by performing an exact diagonalization of the many-body 
Hamiltonian 



31 32 33 34. 



(9) 



3 



The kinetic energy term 



Wi = S + Ej 



(10) 



includes the effects of a neutralizing background. Here 



S 



el \f2i 



i 

E 

fei ,fe2 



(1 - erf(v^)) 



(11) 



is the Madelung energy, e the electron charge, e the di- 
electric constant, z = ir (\ 2 kl + fcf)/A, and the prime 
in the summation means that the term k\ = ki = is 
excluded. The expression for the scattering matrix ele- 
ment Ai 1 i 2 i 3 i i depends on the quantum numbers 2112*324, 
where, again, ik — (sk,jk, &k)- For the case of positively 
polarized "spins" (i.e., <7fc = + for all k) we have: 



-4*1 



12*3*4 Wi= + 



5'; 



1 e 1 



h+h,h+j4 2do V 2 



A, 



OO OO 



^ ^ y/k* + A 2 fc| Vm 



fci— 1 k.2 — — 00 



cos 



—h(ji -33) 



si! 



s 3 ! 



S4 



52 1 / 1 N ltl+*2+pi+P2 

^1 ^ ( Sl ^ ( S 3 - *») ! V *1 + 2 Pl A *2 + 2 P2 ' 

tl=0 pi = l *2=0 P2 = l 



x 2* 1+t2 - 1 k 2 



2vrA 
m 



2(ti+pi+t 2 +P2 




2(P1+P2) 



A? 



(12) 



where 

£S(x) = £(-ir( n+a 

n \ 1 \ ) 1 n _ m ; j 

m=0 v 7 

are the Laguerre polynomials and 

c _/(s4-*i) if(s 4 -*i)<0 
UPl \ Int{(s 4 - £i)/2} otherwise 

and the same yields for up?,, with 54 — > S2 and ii — » £2, 
[Int{a;} is the integer part of a;]. 

Due to the presence of the spinors cr, the two branches 
ip and the coupling of the Landau levels in pairs, 
the derivation of the complete expression for Ai 1 i 2 i 3 i i 
is highly nontrivial. Moreover, the Hamiltonian matrix 
can be very large and its diagonalization computationally 
very intensive. We calculated the ground state energy per 
particle Eq for a system with four electrons in the lowest 
two Landau levels with different filling factors v = 4/iV s , 
where N s = 8, 9, . . . , 20. The results are shown in Fig. 
for B = 10 Tesla. We see that the presence of SO cou- 
pling lowers considerably the value of Eq, compared to 
the result for a = 0, which coincides with the usual re- 
sults obtained for the FQHE with no SO, with both the 
Zeeman and the kinetic energies included. 

The most intriguing effect of SO coupling in 2DEG, 
however, is found in connection with the magnitude of 
the quasiparticle-quasihole energy gap E gi derived from 
the positive discontinuity of the chemical potential at the 
filling factor v 0. 



As shown in Fig. [5] for four electrons and a filling fac- 
tor of i a nd in Fig. |3| for even smaller values of v = jr 
and j |l0j . large values of a cause the enhancement of 
Eg. This enhancement is larger for small magnetic fields 
(B ~ 1 Tesla, i.e., fields for which the Rashba term is 
still comparable to the Zeeman term), and can be of the 
order of 25% for v = and mid-range values of the cou- 
pling constant (a — 20) . This is seen in Fig. [21 where E g 
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FIG. 2: Quasiparticle-quasihole energy gap as a function of 
magnetic field B for a filling factor v — 1/3 and three different 
values of the SO coupling constant a — 0, 20, 40 calculated for 
four electrons in the lowest two Landau levels in InAs. 
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increases from 0.20 meV for a — to ~ 0.25 for a = 20. 
Smaller increases occur at smaller filling factors, where 



B (Tesla) 



FIG. 3: Quasiparticle-quasihole energy gap as a function of 
magnetic field B for a filling factor v = 1/5 and 1/7 and 
two different values of the SO coupling constant a — 0, 40 
calculated for three electrons in the lowest two Landau levels 
in InAs. 



Eg shows a 17% and 7% increase for v 



T5 and A"' 



respectively. Reasons for this behaviour derive from a 
complex interplay between the different one- [Eq. QlOp] 
and two-body [Eq. (|12[) ] terms. For a ^ 0, one clear 
effect of the SO interaction is that the kinetic energy is 
no longer constant for a given value of the magnetic field 
[Eq. (7)] but depends on the strength of the SO coupling. 
This, coupled with the fact that the interaction terms are 
profoundly modified by the SO interaction, results in a 
change in dE / dv. This is reflected in an increase of the 
energy gap. 



In summary, we have investigated the influence of the 
SO coupling (Bychkov-Rashba) on the incompressible 
state proposed by Laughlin at v = -|, using the exact 
diagonalization scheme for finite-size systems in a peri- 
odic rectangular geometry. We found that, as the SO 
coupling strength is increased, there is an increase of the 
quasiparticle-quasihole gap. This is particularly advan- 
tageous for filling factors v < 1/5, 1/7 where a larger gap 
would signify a more stable liquid state that will push 
the liquid-solid transition further down in the density. 
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